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Abstract 

CDM models with non-scale-free step-like spectra of adiabatic pertur- 
bations produced in a realistic double inflationary model are compared 
with recent observational data. The model contains two additional free 
parameters relatively to the standard CDM model with the flat (n — 1) 
initial spectrum. Results of the COBE experiment are used for the deter- 
mination of a free overall spectrum normalization. Then predictions for 
the galaxy biasing parameter, the variance for "counts in cells" , the galaxy 
angular correlation function, bulk flow peculiar velocities and the Mach 
number test are obtained. Also considered are conditions for galaxy and 
quasar formation. Observational data strongly restricts allowed values for 
the two remaining model parameters. However, a non-empty region for 
them satisfying all considered tests is found. 



1 Introduction 

Inflationary cosmological models (Starobinsky 1980, Guth 1981, Linde 1982, 
1983) imply a density parameter £l t ot ~ 1 (| ^tot — 1 |< 10~ 4 ) within the 
observable part of the Universe. Combined with the value f^ar « 0.017/i~ 2 
following from the theory of primordial nucleosynthesis (see, e.g. Walker etal. 
1991) this requires the most of matter in the Universe to be nonbaryonic (h — 
-ffo/100 km/s/Mpc). Second, the simplest inflationary models (with one slowly 
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rolling effective scalar field) predict density fluctuations with a spectrum of 
approximately Zeldovich-Harrison type (i.e., (y)^ oc k n with n w 1). This 
result was first consistently derived by Hawking 1982, Starobinsky 1982, and 
Guth & Pi 1982. The cold dark matter (CDM) model with this initial spectrum 
of perturbations and with biasing b g ~ 1.5 — 2.5 have most successfully explained 
the observed hierarchy of cosmic structures up to scales of approximately 10 
h^ 1 Mpc (Davies etal. 1985). However, observations of such structures as the 
Great Attractor (Lynden-Bell et al. 1988, and Dressier 1991) and the Great Wall 
(de Lapparent etal. 1986), the large-scale clustering in the redshift survey of 
IRAS galaxies (Efstathiou etal. 1990a, Saunders etal. 1991), and the galaxy 
angular correlation function and "counts in cells" for the deep APM galaxy 
survey (Maddox et al. 1990, see also Loveday et al. 1992) imply that there is 
more power in the perturbation spectrum at scales larger than approximately 
10/i -1 than expected in the standard model (biased CDM plus n = 1 initial 
spectrum of adiabatic perturbations). On the other hand, it is remarkable 
that the standard model is so close to these observational data: only a modest 
increase in amplitude of perturbations at large scales is required to fit them - no 
more than 2-3 times, and the latest results have the tendency to diminish this 
number (see e.g. Loveday etal. 1992 ). Finally, recent COBE measurements 
(Smoot etal. 1992) also imply the ratio of amplitude at scales (10 3 — 10 4 )/i _1 
Mpc to that at (1 - 10)/i _1 Mpc equal to (1.1 ± 0.2)b g (la error bars). The 
earlier published positive RELICT - 1 result for AT/T at large angles (Strukov 
et al. 1992a, b) is even larger but it was obtained using only one wavelength, so 
it may be at least partially non-primordial. 

The spectrum of perturbations observed at the present time is a product of an 
initial (primordial) spectrum and some transfer function T(k). The latter results 
from a transition from the radiation-dominated era to the matter dominated one 
at redshifts z ~ 10 and depends on the structure of dark matter. Thus, any 
deviation from the standard model may be explained either by changing the 
initial flat spectrum, or by having a different T(k) due to a more complicated 
matter content. The latter possibility arises, for example, in models with a 
mixture of hot and cold dark matter (Shafi & Stecker 1984, Holtzmann 1989, 
van Dalen & Schaefer 1992, and others) or with a cosmological constant and dark 
matter (Peebles 1984, Kofman & Starobinsky 1985, Efstathiou etal. 1990b, and 
Gorski etal. 1992). For a recent reanalysis of the former model see Pogosyan 
& Starobinsky 1993; the same for the latter model - see Bahcall etal. 1993. 

Here, we consider CDM models with non-scale-free primordial perturbation 
spectra. The scale invariance of the perturbation spectrum can be broken by 
different physical mechanisms during an inflationary stage (see, e.g., Kofman 
et al. 1985, Salopek et al. , 1989). Scale-free but not scale-invariant spectra with 
n < 1 seem not to be able both to provide enough power at L ~ (25 — 50) ft. -1 
Mpc and to fit the COBE data (Polarski & Starobinsky 1992, Liddle etal. 
1992, Adams etal. 1993; see, however, Cen etal. 1992). A non-scale-free initial 
spectrum with an effective step ( compared to the n = 1 spectrum) somewhere 
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between l/i -1 Mpc and lOh^ 1 Mpc works much better. 

To obtain such a spectrum, it is necessary to abandon at least one of the 
main assumptions of the simplest version of the inflationary scenario, i.e. to 
assume either that the slow-rolling condition is temporarily violated at some 
moment of time (and then a step of an universal form in the spectrum arises, 
see Starobinsky 1992), or that there is more than one effective scalar field driving 
inflation. The latter possibility leads to double inflation , i.e., to a cosmological 
model with two subsequent inflationary stages (Kofman etal. 1985, Silk etal. 
1987). These inflationary stages may be driven by the R 2 term and a scalar 
field (Gottlober ct al. 1991), where R is the Ricci scalar or by two noninteracting 
scalar fields (Polarski & Starobinsky 1992). In this case, a step in the spectrum 
arises as a result of a rapid decrease of the Hubble parameter H = a/ a in a 
period between two inflationary stages with slowly varying H and H. Further, 
we consider the first of the above mentioned double inflationary models, results 
for the second one are expected to be similar. 



2 Perturbation spectrum and normalization 

The Lagrangian density of the gravitational field including the R 2 term and a 
massive scalar field reads 

The R 2 term is coupled via M 2 <C G _1 to General Relativity, (p is the scalar 
field and m 2 <C G _1 is its mass (c = fi = 1). The model contains two more 
parameters than the standard model. They define the location and the relative 
height of an effective step in the perturbation spectrum. The height of the step 
depends on the ratio M/m and the energy density of the scalar field at the onset 
of the second inflation. On the other hand, this energy density is responsible for 
the length of the second inflationary stage, i.e., for the physical scale at which 
a break in the spectrum appears. With if = tp ~ 3G~ 1//2 at the end of the first 
inflationary stage, the break occurs just at the right place, and the height of 
the step is A « M/6.5m. Note that due to the exponential dependence of the 
break location on the energy density, its shifting from 10 Mpc to 100 Mpc does 
not practically influence the height of the step. 

We consider the rms value of the Fourier transform of an gravitational 
perturbation <I>(fc). Here the conventions <I>(fc) = (27r)~ 3 / 2 / §(f)e~ lkif d 3 r, 
($(£)) = 0, ($(fc)$*(fc')> = $ 2 (k)6 3 (k - k'), k = \k\ for initially Gaussian per- 
turbations are adopted. Provided the limit k — > of <P(k) coincides with a flat 
spectrum <P(k). Then we can construct the two quantities = lim^fc^-o <P(k) 
and <P 2 = lim fe _ 

febr+o &{k). The height of the step A is given by A= <pW/$( 2 ). 
In this paper we consider spectra with steps A w M/6.5m = 2, 3, 4, 5. The 
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Figure 1: Density perturbation spectra Pj (i =1,2,3,4) for M = 13m, M = 
20m, M — 26m, M = 32.3m respectively, normalized to (7t(10 o ) by the use of 
the COBE data. The dashed line shows the standard CDM model with the flat 
initial spectrum normalized in the same way. 

step is located at L ~ 2-Kk^}, where k^ 1 denotes the wave number where the 
perturbation spectrum reaches the lower plateau at large k, the values k^} = 
(1,3,7,10,20,30) Mpc -1 were considered. Here and further throughout the pa- 
per, unless a scaling by h is given, we assume h — 0.5. The length of a transition 
period between flat parts of the spectra is of the order of one magnitude of wave 
numbers. 

Results of numerical calculations of the power spectra of density perturba- 
tions at the present epoch 

P(k) ee P(k, z = 0) = (^j = l(fci? ff ) 4 $ 2 (fc)T 2 (fc) (2) 

are plotted in Fig.|l|. Here T(k) is the transfer function for the standard CDM 
model (a review of best fits to it see in Liddlle & Lyth 1993), for equations for 
perturbations see Gottlober etal. 1991, the overall normalization is explained 
below. We adopt Sltot = 1 in accordance with inflation, thus a(t) oc t 2 / 3, now. 
Rh = ZHq 1 ~ 6000/i _1 Mpc is the present day Hubble radius sometimes non- 
rigorously called the cosmological horizon. For A < 3, there is practically no 
oscillations in the spectrum. Then a very good analytical approximation for 
it may be obtained by assuming that the slow-roll condition is valid for both 
inflaton fields R and 4> in the region of the break (see, e.g., Gottlober etal. 
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1991). 

Let us present this solution in a different form. The background equations 
for the model (1) in the slow-roll approximation 



have an integral 



3H(f> + m 2 (j) = 0, (3) 
~W = 3 H (4) 



a 3H 2 



ln-- = 2^ + ^_, (5) 

where a/ is the value of the scale factor at the end of inflation (Kofman etal. 
1985, Starobinsky 1985a). Then, introducing / = H 2 /s, we arrive to the equa- 
tion: 



'" =-Af-.h){f2-.f), (6) 



dins / 

a = \™\ h = \m\ m < h- 

The scalar field follows from the relation: 

2irG<p 2 =s(l- ^ . (7) 
Let, e.g., M > my/2. Then the solution of Eq.(6) is 

{f2 - f)7 T =>-/!)• W 

f may be both larger and smaller than f\. f — for s = 2s\. For values of 

the quantities M/m and used in our calculations, f\ <C ji and s\ » 1. Using 
the general expression for adiabatic perturbations generated during a slow-roll 
multiple inflation (Starobinsky 1985a) we get the answer: 

fc3$ 2 (fc) = JL . l67rGs 2 /(l - -4"). (9) 
w 100 Jy 2f 2 

where s and f(s) are taken at the moment of horizon crossing at the inflationary 
stage k = a(t)H(t); so s — s(k) w In in this expression, fc/ = ajHf. 

It is useful to note here that the background dynamics of an inflationary 
model with two massive scalar fields (see Polarsky & Starobinsky 1992) in the 
slow-rolling regime may be reduced to the same equation (6) with / = H 2 / s 
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and /1.2 = f m i 2 ( m i < TO 2)- Thus, the replacement rule is A/ 2 — > 2m 2 ,. 
The integral (7) is then changes to s = 2ttG(4>\ + 4>\). fx < f < fi in this 
case, so inflationary dynamics of the model with two massive scalar fields is 
isomorphous to only a finite part of possible inflationary regimes of the model 
(1). The expression for the power spectrum of adiabatic perturbations is slightly 
different from (9): 

fc 3 $ 2 (fc) - ^ • 16ttG S 2 /. (10) 

Now let us turn to an overall normalization of the power spectra. It is chosen 
to fit the COBE data (Smoot etal. 1992) on large-angle AT/T anisotropy 
(2 < I < 30, where I is the multipole number). If a; is a rms multipole value 
summed over all m and averaged over the sky: 

then the expected variance of the COBE data (Tj,(6fwhm) may be expressed 
as 

°t(0 fwhm ) = ex P + !) 02 ) . (!2) 

Z>2 

where # s = 9fwhm /2\/ln 4 is the Gaussian angle characterizing smearing due 
to a finite antenna beam size, as well as due to an additional Gaussian smearing 
of raw data. We take cr T (10°) = (30 ± 7.5)(j,K/2.735K (Smoot etal. 1992), so 
6 a = 4.25° w 1/13.5. 

We have numerically calculated the rms multipole values a; for the pertur- 
bation spectra plotted in Fig.[| (Gottlober&Mucket 1993). In addition, it is 
possible to obtain an analytical expression for a/. Indeed, if f\ -C / <C then 
Eq.(8) simplifies to / w f 2 (s — Si)/si for < s — si <C si- Thus, si f» ln^-. 
Then it follows from Eq.(9) that 

fc3$ 2 (fc) _9_ . 16 G/ / _ ) = _9_A 2 In ^ (13) 
V ; 100 J u ' 100 k { ' 

for kb r e~ Sl <C fc <C fcbr ( here we introduce the quantity A used in Starobinsky 

1983). Using the standard formula for the Sachs- Wolfe effect which is the only 

significant one for 2 < / < 30 and the case of adiabatic perturbations, it is 

straightforward to derive that 



A 2 (2l + 1) r j2 ri r> Mhr/k) ^ 

A 2 21 + 1 (, , l + \ 



In(2l ) -!-*(/)- —±- , (14) 



400tt 2 Z(Z + 1) V J(' + 1) 



G 



where lo — kbrRn 3> 1 is the multipole number corresponding to the step 
location (/ <C lo) and ^> is the logarithmic derivative of the r-function, the 
difference between Rh and the radius of the last scattering surface may be 
neglected here. Substituting af into Eq.(12) and choosing A to get the correct 
value for or (10°), we obtain properly normalized power spectra P(k) which are 
plotted in Fig.|l]. 

It should be noted that, for the model considered, the contribution of pri- 
mordial gravitational waves (GW) to large-angle AT/T fluctuations is small as 
compared to the contribution from adiabatic perturbations (AP) and may be 
neglected. Really, using formulas derived in Starobinsky 1985a, b, one gets 



a J^l „ 2 .5 ./— * 0.1 (15) 

for 1 <C I < 30, the value of this ratio for I = 2 being « 6% more, fc# = 2ttR^ h 1 . 
Thus, the account of the GW contribution increases the total rms value of 

ai(tot) = J a2 i(AP) +a2 i(GW) b y less than 1% - 



3 Comparison with observational data 
3.1 Biasing parameter 

By help of the density perturbation spectrum P(k) we are able to compute the 

variance of the mass fluctuation <j 2 M (R) = (-^f(R)) 2 in a sphere of a given 
radius R (see e.g. Peebles 1980): 

1 f°° 

<J 2 M (R) = ^J q k 2 P(k)W(kR)dk, (16) 
where the window function W(kR) is expressed as 

W(kR) = -j^isinkR- kRcoskR) 2 (17) 

The quantity a 2 M gives the variance of total matter density perturbations. It 
becomes clear now that the distribution of galaxies is biased with respect to 
the distribution of total matter, so some biasing parameter b g (assumed to be 
scale-independent for simplicity) has to be introduced: 

Ur) = b 2 g m°M g (R) = by M (R), (18) 

where £(r) is the matter density correlation function. Usually it is adopted 
that cr1 Ig (8h^ 1 Mpc) = 1, so b g = a^l^h^ 1 Mpc). Comparison of N-body 
simulations with the observed £ g (r) and with the mean pairwise velocity of 
galaxies at r = lh^ 1 Mpc (Davies etal. 1985, for more recent calculations 



7 



Spectrum Pi[k^ r L ] 


Bias factor 


(Tv(R) 


(Tv(R) 


depending on step 


b 9 


km/s 


km/s 


location at kz} Mpc 

or L 




R = 40/1" 1 Mpc 


R=60h' 1 Mpc 


Pl(03) 


1.60 


269 


221 


P2(03) 


2.18 


257 


213 


P3(03) 


2.43 


259 


215 


P4(03) 


2.19 


268 


221 


Pl(07) 


1.74 


255 


211 


P2(07) 


2.77 


230 


195 


P3(07) 


3.49 


229 


195 


P4(07) 


3.68 


242 


204 


Pl(10) 


1.77 


248 


205 


P2(10) 


2.83 


215 


184 


P3(10) 


3.59 


210 


182 


P4(10) 


4.18 


225 


193 


Pl(20) 


1.76 


234 


195 


P2(20) 


2.67 


182 


159 


P3(20) 


3.46 


168 


151 


P4(20) 


4.48 


180 


162 



TABLE 1: The bias factors and the velocity variances for the models. PI, P2, 
P3 and P4 denote models with M = 13m; 20m; 26m; 32.2m corresponding to 
A =2, 3, 4, 5. The step location is given by k^ 1 . 

see e.g. Gelb 1993) shows that b g probably lies between 2 and 2.5 (0.4 < 
<tm(8/i _1 ) < 0.5). We adopt rather conservative limits 1.5 < b g < 3 (0.33 < 
c7m(8/i _1 ) < 0.67). Different bias factors obtained for the considered spectra 
are given in Table 1 (one should not forget about 25% error bars at 1.5a level 
in all these quantities due to the error bars in ctt(10 )). This test excludes both 
the standard model and the P(4) case (apart from the region fc^ 1 w 3 Mpc for 
the latter). The P(3) case is marginally admissible, and a broad range of k\, r 
remain allowed for the P(l) and P(2) cases. 



3.2 Counts-in-cells 

Further we compare the variance o\(X) for the counts-in-cells analysis of large- 
scale clustering in the redshift survey of IRAS galaxies (Efstathiou et al. 1990) 
and in the Stromlo-APM survey (Loveday etal. 1992) with our predictions 
obtained for corresponding scales. The variance a^(l) is related to the two- 
point galaxy correlation function according to 
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Figure 2: The variance a 2 (l) plotted as a function of cell size I. Points obtained 
from the counts-in-cells analysis (Loveday et al. 1992) are plotted with 2a error 
bars. Lines show predictions of a 2 (I) computed from our spectra with the step 
location at k^} = 7 Mpc. The solid lines correspond to spectra with different 
step height as indicated at the lines, the dashed line corresponds to the standard 
("flat") CDM model. The IRAS data are plotted as quadratic symbols. 



°2(0 = ^2 J j y _ 3 MrvfidVim . (19) 
In terms of the density perturbation spectrum P(k), it can be written as 

1 r°° 

°^(0 = 2^y vpiQw^tydk, (20) 

where 

W,(p)=8 fdx fdy f dz(l - - y)(l - z ^^ff+f +^l. (21) 
Jo Jo Jo p\f x z + y z + z z 

A comparison between a 2 (I) for galaxies from the Stromlo-APM survey and 
IRAS galaxies (see Fig.||) shows that the biasing parameter for IRAS galaxies 
biRAS < bjipM for I = 10/i -1 Mpc and bmAS ~ bAPM for larger scales (apart 
from the point I = 40/i _1 Mpc where b 2 IRAS w l.bb 2 APM that is, however, inside 
2cr error bars). Predictions for a 2 {l) for the considered spectra normalized by 
the condition (7^j g (8h^ 1 Mpc)=l are shown in Figj^. It is seen that another 
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Figure 3: \ 2 test of variances. The contours in the fc^-A-plane represent the 
levels x 2 (A,/c b ; 1 ) = 3, 7, 20, respectively. 

normalization condition cr c (12.5/i _1 Mpc)=l proposed in Pogosyan & Starobin- 
sky 1993 works as well. The Kaiser correction (Kaiser 1987) was not taken into 
account. In any case, it is clear that this correction should be rather small that 
presents one more argument for a sufficicnly large value of b g . 

In order to determine values of the model parameters which give the best fit 
to these data, we formally apply the x 2 test considering the error bars given in 
Loveday etal. 1992 as 2a ones (in the logarithmic scale of cr c (l)). The results 
are plotted in Fig.|^. Taking into account the possibility of a change in the 
overall normalization, the estimated number of degrees of freedom is N = 7, so 
we consider fits with \ 2 < 7 as good and fits with x 2 > 20 as unacceptable. 
The best fit X 2 ~ 2.2 occurs for A = 3 and k^ 1 = 3 Mpc. 

3.3 Angular correlation function 

Let us consider the predictions of our model for the galaxy angular correlation 
function w(6). To this aim we put the Fourier transform of the 3-D correlation 
function £ g (r) into the Limber equation and find 

, a s Jo°° kP(k)dkJ 00 y^(y)Jo(key)dy 

27T 2 (j y 2 ny)dyy 

We insert our power spectra into Eq.(22) and scale the angular correlation 
function to the Lick depth (y* ss 240h~ 1 Mpc) with the selection function 
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Spectrum Pi\k h ,„ L ] 

" h or J 


Mach numbers 


depending on step 


R L = 4/i~ i Mpc 


R L = 8/i _i Mpc 


R L = 18/i _i Mpc 


location at kz} Mpc 

or * 








Pl(03) 


2.76 


1.48 


0.76 


P2(03) 


3.16 


1.68 
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1.76 


0.89 
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3.21 
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0.85 


Pl(07) 


2.80 


1.51 


0.79 


P2(07) 


3.45 


1.88 


0.98 


P3(07) 


4.01 


2.17 


1.10 


P4(07) 


4.18 


2.19 


1.07 


Pl(10) 


2.78 


1.51 


0.79 


P2(10) 


3.37 


1.87 


1.01 


P3(10) 


3.90 


2.20 


1.19 


P4(10) 


4.51 


2.44 


1.22 


Pl(20) 


2.68 


1.46 


0.77 


P2(20) 


2.85 


1.61 


0.93 


P3(20) 


3.09 


1.76 


1.07 


P4(20) 


3.87 


2.25 


1.35 



TABLE 2: The Mach number test (notation as in Table 1) 



4> = y~ - 5 exp (— (y/y*) 2 ) ■ Results of the calculations along with the recent 
observational data for w(9) for the APM survey (provided by S. Maddox, see 
also Maddox etal. 1990) are shown in Fig.||. The use of linear approximation 
in the calculation of w(9) is justified for 9 > 0.5° only. This is the reason of a 
discrepancy between our theoretical curves and the data at small angles that 
should not be taken into account. At larger angles up to 9 ~ 10°, the agreement 
with the A = 3 and A = 4 cases is very good (and errors are too large to make 
any definite conclusion beyond this angle). 

3.4 Large-scale bulk flows 

Next we investigate our spectra with respect to large-scale bulk flows (peculiar 
velocities) of matter. The rms peculiar velocity inside a sphere of a radius R is 
given by the expression: 

2 f°° 

a li R ) = ^m- P(k)W(kR)exp{-k 2 R 2 s )dk, (23) 

77 R H JO 

where W(kR) is the window function given in eq. ( |l7j ) and the length scale Rs 
characterizes the Gaussian smoothing of raw observational data. In the Table 1, 
the corresponding predictions for our model are shown. They may be compared 
with recent data from the POTENT reconstruction of the 3-D peculiar velocity 
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field: o-„(40/j- x Mpc) = (405±60) km/s and a^Oh^Mpc) = (340±50) km/s 
with Rs = 12/i _1 Mpc, error bars are la (Dekel 1992, see also Bertschinger 
etal. 1990). 

This test presents the most serious problem for the model because expected 
values of a v are lower than in models having exactly flat spectrum at large 
scales, e.g. in the CDM+HDM n = 1 model or in the CDM model with a step 
in the initial spectrum produced by one scalar field (Starobinsky 1992), due 
to the logarithmic increase of P(k) at k <C fc{, r , see Eq. (13). On the other 
hand, for k^ 1 < 7 Mpc, they are higher than in tilted models with n < 0.75 
and without significant contribution of gravitational waves to large-angle AT/T 
fluctuations. It should be mentioned here that bulk flow velocities obtained from 
observations contain large uncertainties. Second, it is fairly possible that actual 
values of bulk flow velocities at our location in the Universe are larger than 
average ones ("cosmic variance"). So, taking the data at their lower 2a limit 
(this is ~ 1.4 times less than the average values that is permitted by cosmic 
variance) and pushing the overall normalization of the model to the upper 1.5a 
limit of the COBE data, we find that the allowed region is k^ 1 < 7 Mpc. 



3.5 Mach number test 

Now we consider the Mach number test for our model (Ostriker & Suto 1990). 
The expression for the Mach number in terms of the quantities introduced above 
is given as 



M 2 (R s = Jo°° m ex P (-fc 2 i?|) eM-k 2 R 2 L )dk 

1 L> J oo P(fc)exp(-fc 2 i?|)(l-(l + A: 2 i?2/9)exp(-fc 2 i?2))dfc' 1 ' 

where Rs = 5/i _1 Mpc characterizes the Gaussian smoothing. Ostriker & Suto 
(1990) found the Mach numbers 4.2 ± 1.0, 2.2 ± 0.5 and 1.3 ± 0.4 at the scales 
Rl = 4/i _1 Mpc, 8/i _1 Mpc, 18/i _1 Mpc, correspondingly. The Mach numbers 
predicted for our spectra are given in Table 2. Here, once more, the values of 
M obtained from observations in our vicinity may be larger than average ones 
due to the numerator of Eq. (|2~4|), i.e., due to local values of large-scale bulk 
velocities being larger than average. This test clearly shows that A > 2. 



3.6 Quasar and galaxy formation 

Finally, we consider the compatibility of our model with the existence of a suffi- 
cient number of large galaxies at the redshift z = 1 and host galaxies for quasars 
at z = 4. The standard model of quasar formation assumes that they arise as 
a result of formation of massive black holes (M sa 10 9 M Q ) in nuclei of galaxies 
with total masses (10 11 — 1O 12 )M (incuding the dark matter component). We 
estimate /(> M) - the fraction of matter in gravitationally bound objects with 
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Spectrum Pi\k h J~] depending on step 


x 2 


location at kZs Mpc 

or L 
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18.00 


Pl(03) 


13.41 


P2(03) 
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P3(03) 


11.08 


P4(03) 


10.17 


Pl(07) 


14.44 


P2(07) 


11.18 


P3(07) 


13.69 


P4(07) 


28.06 


Pl(10) 


15.80 


P2(10) 


14.27 


P3(10) 


13.69 


P4(10) 


20.04 


Pl(20) 


19.56 


P2(20) 


49.82 


P3(20) 


24.17 


P4(20) 


17.27 



TABLE 3: x 2 -values for the models 



masses beginning from M and higher - using a very simple though rather crude 
formula by Press & Schechter (1974): 



'fe^-'-^WK^iJ- (25) 

where M = ^irR 3 po, <tm(R,z) — CTAf(i?)/(l + z), <tm(R) is the rms value of 
the mass fluctuation in the linear approximation at the present moment, it is 
dcfmcd in Eqs. (16,17). The choice of the quantity 5 C that is equal to the value 
of ^ in the linear approximation at the moment when the considered region 
collapses in the course of a fully non-linear evolution is critical for this approach. 
For the spherical model, S c — 3(127r) 2 / 3 /20 ~ 1.69, but attempts to fit results 
of N-body simulations to the Press-Schechter formula did not produce a unique 
answer. It seems that the best fit for S c lies between 1.33 and 2. 

If we accept the most recent estimates of the mass fraction in hosts of quasars 
at z = 4 (see, e.g. Haehnelt 1993): /(> 1O U M ) « 10" 4 , when it follows 
from Eq. ( p5| ) that the corresponding rms linear mass fluctuation at the present 
moment 

CT M (10 n M©) 1.35 c w 2.2 ± 0.5 (26) 

for S c lying in the range mentioned above. Clearly, this is only a lower limit on 
(7m because quasars may form not in all host galaxies and the presenly observed 
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Figure 4: Predicted angular correlation functions compared with the observa- 
tional data from the APM survey (the data were kindly provided by S. Maddox). 
The solid lines correspond to a spectrum with = 7 Mpc, the dashed line 
shows the prediction for the standard CDM model 

quasar density at z = 4 may be less than the real one. However, due to extreme 
sensitivity of the expression ( p5| ) to a change in <tm, the actual value of the latter 
cannot be significantly larger than that given in Eq. ( |26| ) . A similar estimate was 
presented in recent paper by Blanchard etal. (1993) (it is 30% higher for their 
value a = 1.7). 

Another estimate may be obtained from the fact that large galaxies (or, at 
least, a significant part of them) seem to be already existent at z — 1. Then, 
assuming /(> 1O 12 M ) > 0.1 at z = 1, we get 

cr M (10 12 Af Q ) > 1. 28 c w 2.0 ±0.4. (27) 

Note also that the estimate of Haehnelt (1993): /(> 1O 12 M ) > 10~ 5 at z = 4 
- leads to almost the same result: t7M(10 12 M Q ) > 1.1S C . 

Now the same quantities may be calculated for our spectra. Analysis of the 
results shows that the predicted values of <jm (10 n Mo) and o-m(1O 12 M ) for 
A = 3 and 1 Mpc < < 20 Mpc lie just inside the required ranges. For 
A = 2 they are probably too high, for A = 4 - too low (apart from the points 
k^ 1 = 1 Mpc and 20 Mpc). 
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4 Conclusions 



We have compared the double inflationary model (1) characterized by two ad- 
ditional free parameters as compared to the standard model with a number of 
observational tests and have found that there is a reasonable agreement with all 
considered tests for the following region of the parameters: 1 Mpc < fc^ 1 < 10 
Mpc (that corresponds to real length scales Lb r w (6 — 60) Mpc) and 2 < A < 4. 
The best fit seems to be given by = (3-7) Mpc, A 3. Generally, there is a 
correlation between a step size and a break location giving a good fit to the data, 
i.e., for a break at larger scales also a higher step is necessary. However, due 
to a complicated structure of the obtained spectra (e.g., oscillations throughout 
the range of the break for A > 2) this connection is not so straightforward. 

The most restricting for the model are large scale bulk flows. To obtain even 
a marginal agreement we have to assume that their actual rms values are at 
least w 1.4 times less than the average values of Dekel's (1992) data (due to 
the cosmic variance) and that the rms value of large-angle CMB temperature 
fluctuations is not less than the average result of COBE. These assumptions are 
certainly crucial tests for the considered model. 

The work of one of the authors (A.S.) was supported by the Russian research 
project " Cosmomicrophysics" . A.S. also thanks the Deutsche Forschungsge- 
meinschaft (DFG) for the finacial support of his visit to the Potsdam Astro- 
physical Institute. We would like to thank Steven Maddox for sending us the 
latest version of his w{6) data. 



15 



References 



Adams, F.C., Bond, J.R., Freeze, K., Frieman, J. A., and Olinto, A.V., 1993, 
Phys.Rev D47, 426 

Bahcall, N., Gnedin, N.Yu., and Kofman, L.A., 1993 Ap J 413, 1 
Bertschinger, E., Dekel, A., Faber, S.M., Dressier, A., and Burstein, D., 1990, 
ApJ 364, 370 

Blanchard, A., Buchcrt, T., and Klaffl, R., 1993, A&A,267,1 

Cen, R., Gnedin, N.Yu., Kofman, L.A., and Ostriker, J. P., 1992, ApJ Letters 

399, Lll 

Davies, M., Efstathiou, G., Frenk, C.S., and White, S.D.M., 1985, ApJ 292, 
371 

Dekel, A., 1992, in Observational Cosmology, ed. by G. Chincarini etal. , 
Milan 

de Lapparent, V., Geller, M., and Huchra, J., 1986, ApJ Letters 202, LI 
Dressier, A., 1991, Nature 350, 391 

Efstathiou, G., Kaiser, N., Saunders, W., Lawrence, A., Rowan-Robinson, 
M., Ellis, R.S., and Frenk, C.S., 1990a, MNRAS 247, lOp 
Efstathiou, G., Sutherland, W.J., and Maddox, S. J., 1990b, Nature 348, 705 
Gelb, J.M., 1993, ApJ , in press 

Gorski, K.M., Silk, J., and Vittorio, N., 1992, Phys. Rev. Lett. 68, 733 
Gottlober, S., Miicket, J.P., 1993, A&A 272, 1 

Gottlober, S., Mullcr, V., and Starobinsky, A.A., 1991, Phys. Rev. D 43, 
2510 

Guth, A.H., 1981, Phys. Rev. 23, 347 

Guth, A.H., and Pi, S.Y., 1982, Phys. Rev. Lett. 49, 1110 

Hawking, S.W., 1982, Phys. Lett. B115, 295 

Haehnelt, M.G., 1993, MNRAS, in press 

Holtzman, J. A., 1989, ApJ Suppl.Series 71, 1 

Kaiser, N., 1987, MNRAS 227, 1 

Kofman, L.A., and Starobinsky, A. A., 1985, Sov.Astron.Lett. 11, 271. 
Kofman, L.A., Linde, A.D., and Starobinsky, A.A., 1985, Phys. Lett. B157, 
361 

Linde, A.D., 1982, Phys. Lett. B108, 389 
Linde, A.D., 1983, Phys. Lett. B129, 177 
Liddle, A.R., and Lyth, D.H., 1993, Phys. Rep. 231, 1 

Liddle, A.R., Lyth, D.H., and Sutherland, W., 1992, Phys. Lett. B279, 244 
Loveday, J, Efstathiou G., Peterson B.A., and Maddox, S.J., 1992, ApJ Let- 
ters 400, L43 

Lynden-Bell, D., Faber, S.M, Burstein, D., Davies, R.L., Dressier, A., Ter- 
levich, R.J., and Wegner, G., 1988, ApJ Letters 326, L9 
Maddox, S.J., Efstathiou, G., Sutherland, W.J., and Loveday, J., 1990, MN- 
RAS 242, 43p 



16 



Ostriker, J.P., and Suto, Y., 1990, Ap J 348, 378 

Peebles, P.J.E., 1980, The Large-Scale Structure of the Universe (Princeton, 
Princeton University Press) 
Peebles, P.J.E., 1984, ApJ 284, 439 

Pogosyan, D.Yu., and Starobinsky, A. A., 1993, MNRAS, in press 
Polarski, D., and Starobinsky A. A., 1992, Nucl. Phys B385, 623 
Press, W.H., and Schcchtcr, P., 1974, ApJ 187, 425 

Salopek, D.S., Bond, J.R., and Bardeen, J.M., 1982, Phys. Rev. D40, 1753 

Saunders, W., et al. , 1991, Nature 349, 32 

Shaft, Q., and Stecker, F.W., 1984, Phys. Lett. 53, 1292 

Silk, J., and Turner, M.S., 1987, Phys. Rev. D35, 419 

Smoot G.F., Bennett C.L., Kogut A., Wright E.L. etal. , 1992, ApJ Letters 
396, LI 

Starobinsky, A.A., 1980, Phys. Lett. B91, 99 
Starobinsky, A.A., 1982, Phys. Lett. B117, 175 
Starobinsky, A.A., 1983, Sov.Astron.Lett. 9, 302 
Starobinsky, A. A., 1985a, JETP Lett. 42, 152 
Starobinsky, A. A., 1985b, Sov.Astron.Lett. 11, 133 
Starobinsky, A.A., 1992, JETP Lett. 55, 489 

Strukov, LA., Brukhanov, A. A., Skulachev, D.P. and Sazhin, M.V., 1992a, 
Sov.Astron.Lett. 18, 153 

Strukov, LA., Brukhanov, A. A., Skulachev, D.P., and Sazhin, M.V., 1992b, 
MNRAS 258, 37 

van Dalen, A., and Schaefer, R.K., 1992, ,4pJ398, 33 
Walker, T.P. et al., 1991, ApJ 376, 51 



17 



